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| Abstract 

£NJ ■ Both classical and quantum damped systems give rise to complex spectra and 

corresponding resonant states. We investigate how resonant states, which do not 
D ■ belong to the Hilbert space, fit the phase space formulation of quantum mechanics. 

^0 \ It turns out that one may construct out of a pair of resonant states an analog of 

the stationary Wigner function. 

> ■ 1. Introduction 

00 

Recently, it was observed by Kossakowski [1], [2] that quantized simple damped sys- 
Q\ ■ terns (e.g. damped harmonic oscillator) give rise to discrete complex spectra. The 

corresponding eigenvectors may be interpreted as resonant states. Such states play im- 
portant role in quantum mechanics and it is widely believed that they are responsible 
for the irreversible dynamics of physical systems (see e.g. [3]). In a recent paper [4] it 
was shown that the damping behaviour in a classical system may be also interpreted 
as appearance of resonant states for the corresponding Koopman operator. Consider a 
classical Hamiltonian system defined by the Hamilton function H : P — > R. Assuming 
that a Hamiltonian flow is complete one introduces a self-adjoint operator Lh in the 
Hilbert space L 2 {P,dp), where d[i denotes the standard Liouville measure on P. This 
^ ■ so called Koopman operator is defined by 

"X' 

where { , } denotes the Poisson bracket in the algebra of classical observables C°°(P). 
Studying simple examples of classical damped systems we showed [4] that corresponding 
Koopman operators have discrete complex spectra. Clearly, the generalized eigenvec- 
tors do not belong to the Hilbert space L 2 (P,d[i) and the appropriate mathematical 
language to deal with this situation is the rigged Hilbert space or the Gelfand triplet 
(cf. e.g. [5], [6]), that is, any generalized eigenvector belongs to a dual space D* where 
D is an appropriate space of test functions in L 2 (P,dfi). 

In the present paper we shall study the phase space formulation of quantum damped 
systems. Clearly, this formulation, called also the deformation quantization, is perfectly 
equivalent to the standard Hilbert space approach. However, as we already mentioned, 
resonant states lie outside the Hilbert space, and hence, it would be interesting to find 
how they fit phase space approach. Any vector tp £ 7i gives rise to a Wigner function 
W^p on a classical phase space P. As was shown already by Wigner [7] this function is 
real and produces marginal probability densities f W^(x,p)dx and J W^(x,p)dp. The 
classical limit of reproduces a classical probability distribution on P. Moreover, if 
ip is an eigenvalue of the Hamilton operator H, i.e. Hijj = Etp, then the corresponding 
Wigner function satisfies the following eigenvalue problem: 

H*Wj, = W,p*H = EWtP , (1.2) 
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L H f:=i{f,H}, (1.1) 
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where H is a classical Hamiltonian on P and / * g denotes "quantum deformation" of 
a usual commutative product of functions / • g (see the next section) . 
It turns out that resonant states appear always in pairs: 

= E ± i/j ± , (1.3) 

such that 

E~ = E+ . (1.4) 
We show that each pair gives rise to a pair of function on P satisfying 

H*F ± = F ± *H = E ± F ± , (1.5) 

with 

F~ = F+ and E~ =F+ . (1.6) 

Moreover, if ^ + = ip~ = 4> G H, then E + = E~ G R and F + = F~~ = , 
that is, F ± may be considered as a generalization of Wigner function for resonant 
states. Functions F^ do indeed satisfy basic properties of Wigner function: they 
may be normalized so that J p F^d^ = 1, and they give rise to marginal probability 
distributions: if (x±, . . . ,x n ,pi, . . . ,p n ) are canonical coordinates on P, then 



Tr^(xi,...,x n ) - J F ± (x 1 ,...,x n ,p 1 ,...,p n )dp 1 ...dp n , (1.7) 

Kp(Pl,---,Pn) ■= J F ± (x 1 , . . . ,X n ,pi, . . . ,p n ) dX! . . . dx n , (1.8) 

are classical probability distributions on P. Actually, in the examples considered in 
this paper one finds 

(xi, ...,x n ) = S(xi) . ..6(x n ) , (1.9) 

ir±(p 1 ,...,p n )=6(p 1 )...5(p n ) . (1.10) 



It seems that the above result does contradict Heisenberg uncertainty relations. Note, 
however, that resonant states do not belong to the Hilbert space and hence the proba- 
bilistic interpretation is not clear. 

In the next section we recall phase space approach to quantum mechanics. Section 3 
shows how this approach works for the harmonic oscillator. Following sections discuss 
damped systems: a toy model of a damped motion x = —jx, and damped harmonic 
oscillator. We end up with some conclusions. 



2. Phase space formulation of quantum mechanics 



Deformation quantization consists in replacing a commutative algebra of functions 
C°°{P) over a classical phase space P by a noncommutative algebra (C°°(P),*). For 
simplicity let us assume that P = M. 2N . The ^-product operation 



★ : C°°(P) x C°°(P) 

is defined by: 

/ * 9 ■= f exp 
where A denotes a bidifferential operator 

f*kg:={f,g} 



— A 
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C°°(P) 



(2.1) 



(2.2) 
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The above structure was introduced long ago by Groenewold [8] and later on it was used 
by Moyal [9] to construct the so called phase-space formulation of quantum mechanics 
(see e.g. recent review by Zachos [10]). The equivalence of the above approach to the 
standard Hilbert space one is based on the well known Weyl correspondence: if A is 
self-adjoint operator on H = L 2 (M. N , dx), then one defines the symbol A(u,v) of the 
operator A by: 

A(u,v) : =Tr(le ( ~ i/?l)(up+vX) ) , (2.3) 

where x = (xi, . . . , xjv) an d p = (pi, . . . ,Pn) are standard position and momentum 
operators in L 2 (R^,e£x). Conversely, given a symbol A(u,\) one construct a corre- 
sponding operator: 

A = J duj dv^(u,v)e(^( u P +vS ) . (2.4) 

Now, j£C = AB, then 

C=A*B, (2.5) 

where A, B, C are symbols of A, B and C, respectively. In this approach the von 
Neumann equation for the density operator p 

ihd t p=[H,p], (2.6) 

is replaced by the Moyal equation for the corresponding Wigner function W, i.e. symbol 
of p: 

ihd t W = {H,W} M , (2.7) 
where the Moyal brackets is given by: 

{H,W} M :=H*W- W*H . (2.8) 
Using (2.3) it is easy to show that Wigner function corresponding to p is given by [7] 

^(x,p):=^ip | dye -w(x-|y|p|x+|y), (2.9) 

where |x) is a generalized eigenvector of x, i.e. x|x) = x|x). This way quantum 
mechanics may be formulated entirely in terms of objects living on a classical phase 
space P. The very definition of the ^-product implies 

f*g-g*f = ih{f,g} + 0(h 2 ) , (2.10) 

and hence in the classical limit h — ► the Moyal equation (2.7) reproduces the Liouville 
equation 

9tU = {H,U}, (2.11) 
for the classical density probability f c \ on P. Now, the unitary evolution of p: 

p(t) = U(t)pU-\t) , (2.12) 

with U(t) = exp(—(i/h)tH), is replaced by the following formula for W(t): 

W(t) = U*{t)*W*U- 1 {t) , (2.13) 

where the so called ^-exponential U± is defined by [12]: 

U^(t) = Exp(-(it/h)H) 

:= 1 + (-it/h)H + ^(-it/h) 2 H * H + ^{-it/hfH * H * H + . . . (.2.14) 
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In the classical limit 

W(t) — / d (t) = e~ itL » /ei , (2.15) 

that is, 

/ c i(x, p, t) = / d (x(-t), p(-t), 0) , (2.16) 

where x(t) and p(t) stand for the classical evolution of x and p, respectively. Actually, 
using (2.1), one finds that the quantum evolution of x and p (in the Heisenberg picture) 

. _ »»H-g** = {xJ/}| (217) 

p = P * g /*" ={p.g}. (2.18) 

is the same as the classical one. 

Finally, let us turn to the energy spectrum. In the standard approach one solves 
for the standard eigenvalue problem for the quantum Hamiltonian H: 

Hip = E^ . (2.19) 

It is easy to see that the corresponding Wigner function W satisfies: 

W$ * H = H * Wj, = EW^ . (2.20) 

Actually, one may prove (see e.g. [13]) that any real solution W of (2.20) corresponds 
to a Wigner function for tp satisfying (2.19). Moreover, if ip n define an orthonormal 
basis in H, then corresponding Wigner functions W n satisfy: 

W n *W m = 5 nm W n , (2.21) 

and hence one obtains the following resolution of identity on P: 

X>-=pSj*- (2 - 22> 

which is phase space analog of the Hilbert space formula J2 n = \ where P n is a 
1-dimensional projector onto the eigenspace generated by ip n . For more properties of 
Wigner function see e.g. the review article [11]. 



3. Harmonic oscillator 

To get the feeling how this approach works in practice let us consider 1-dimensional 
harmonic oscillator described by the following Hamiltonian: 

H ho (x,p) = |( p 2 + x 2 ) . (3.1) 

Now, let us study the corresponding eigenvalue problem (2.20). Equation H*W = EW 
gives: 

2 i 2 
X + p 

whereas W * H = EW: 

2 i 2 

X +p 



2E~ 

UJ 



W = , 



(3.2) 



2E~ 

UJ 



W = . 



(3.3) 



4 



Therefore, 



(xdp-pd x )W = 



(3.4) 



which means that W is a zero-mode of the Koopman operator Lx ho [4]. Taking into 
account (3.2) and (3.3) we obtain: 



* 2 +p 2 -^+3 p 2 )-^ 



W = . 



Introducing a new variable £: 



£:=^ 2 +p 2 ) 



equation (3.5) may be rewritten as follows: 



£ ^2 ~ ^ 



W(0 = . 



Finally, defining L = L(£) by 



W(0 =: e"«/ 2 L(0 , 



equation (3.7) implies: 



£(0 = o 



which is the defining equation of Laguerre's polynomials: 

n! 



(3.5) 
(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 



for n = E/Tluj — 1/2 = 0, 1, ... . This way one recovers well known oscillator spectrum. 
The corresponding Wigner functions W n read 



W n = { -^- e'^ L n (0 . 



(3.11) 



The reader will easily check that W n defined in (3.11) do indeed satisfy formula (2.21). 
It is well known that only Wq which is given by the Gaussian distribution 



1 

irh 



^ = le^ = 4e-( lV '/ fi 



1 

irh 



(3.12) 



defines a probability distribution on P. However, in the classical limit H — ► all 
Wigner functions W n tend to well defined classical probability distributions. For ex- 
ample 

W (x,p) — 6(x)5(p) . (3.13) 

There is an alternative way to find the eigen- Wigner functions W n . One introduces 
phase-space analogs of creation and annihilation operators: 



a = 



x + ip 



x — ip 



f 2h V2h 

satisfying standard commutation relation: 

{a, cl*}m = a* a* — a* *a = 1 



(3.14) 



(3.15) 
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It is easy to rewrite the formula for the ★-product (2.1) in terms of a and a*: 

f- kg = fe 2 \ Jg, (3.16) 

where d = d/da and d* = d/da*. Oscillator Hamiltonian (3.1) takes in the new 
variables the following form: 

H ho = hLo(a**a + l) . (3.17) 



2, 

Now, let us define Wq as a *-Fock vacuum, that is, 

a*W = W * a* = , (3.18) 
and the corresponding excited states: 

W n (x a* n *W *a n , (3.19) 

Noting that £ = 4|a| 2 it is easy to check that W n defined in (3.19) agrees with the 
formula (3.11). 

Finally, let us turn to the time evolution defined in (2.13). The corresponding 
^-exponential (2.14) was found in [12] and is given by: 

Exp(-(z/ft)itf ho ) = 1— exp [-2{i/h) tan(t/2)tf ho ] . (3.20) 

cos(r/2) 

Actually, for the harmonic oscillator the quantum evolution has the same from as the 
classical one. It is evident from (3.2) and (3.3) that 

{H ho ,W} M = ih{H ho ,W} , (3.21) 

and hence, the Moyal equation (2.7) is the same as the Liouville equation (2.11). There- 
fore, due to (2.16) 

W(x,p,t) = W(x(-t),p(-t),0) . (3.22) 

4. Toy model of damped system 

Now, we apply this scheme to the simple damped system described by the following 
equation: 

x = — -yx , (4-1) 

where 7 > is a damping constant. Clearly, this system is not Hamiltonian. However, 
following [14] we may lift an arbitrary dynamics on a configuration space Q 

x = X(x) , (4.2) 

where X is a vector field on Q, to the Hamiltonian dynamics on the corresponding 
phase space P = T*Q. We define the corresponding Hamiltonian 

H : P — ► R , (4.3) 

by 

H(a x ) := a x (X(x)) , (4.4) 
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for a x G T*Q. Using canonical coordinates (x±, . . . , xn,Pi, . . . ,Pn) on P we may 
rewrite a formula for H in a more familiar way: 



A? 



H(x,p) =Y J PkX k (x) . (4.5) 



k=i 

The corresponding Hamilton equations read as follows: 



x k = {x k ,H} = X k (x) , (4.6) 



N 



for A; = 1, . . . , N. In the above formulae { , } denotes the canonical Poisson bracket on 

l".Oi=±(£¥--»°f). (48) 

Clearly, the formulae (4.6) reproduce our initial dynamical system (4.2) on Q. 

Now, applying the above procedure to (4.1) one obtains the Hamiltonian system on 
M 2 with the Hamiltonian given by: 

H d (x,p) = -jxp . (4.9) 

This system was analyzed in [4] where both classical spectrum of the corresponding 
Koopman operator Ln d and quantum spectrum of 

H d = -^(xp + px) , (4.10) 

were found: 

Spec(L^ d ) = {i-fn | n G Z} , (4.11) 

and 

Spec(tf d ) = | ih-y (n + ^ | n € Z | . (4.12) 

Both spectra are discrete and purely imaginary. It should be stressed that both Ln d and 
H d are self-adjoint operators on the corresponding Hilbert spaces L 2 (R 2 ) and L 2 (R), 
respectively. The corresponding eigenvectors, which obviously do not belong to the 
Hilbert space, are usually called resonant states (see e.g. [3]). It was found in [4], [2] 
that for 

^+{x):=x n and ip~(x) := (-ih) n 5 in) (x) , n = 0,l,2,..., (4.13) 

one has: 

H^ = ±ih 7 (n+^J^ . (4.14) 

Evidently, these states living outside the Hilbert space L 2 (M) can not be used to con- 
struct stationary Wigner functions. Indeed, defining 

W£{x,p) oc J dye-^ix - |y) ^{x + \y) , (4.15) 

one obtains 

W£{x,p,t) = e^ 2n+1 ^W,t(x,p,0) , (4.16) 
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which shows that are non-stationary. To find the analogs of stationary Wigner 
functions consider the eigenvalue problem: 

H d *F = EF and F * H d = EF . (4.17) 

Equation H d * F = EF gives: 



X P + -j d x P + y (p 9 p ~ xd x) ~ — 



whereas F * H d = EF: 



n 2 



ih 



E 
7J 



xp + '-^d 2 x - "-^{pd p - xd x ) - — 



Therefore, in analogy to (3.4) one finds 

(pd p - xd x ) F = , 



F = , 



F = . 



(4.18) 



(4.19) 



(4.20) 



which means that F is a zero-mode of the corresponding Koopman operator L}j d . 
Taking into account (4.18) and (4.19) we obtain: 



ft 2 £ 
xp+-d xp + -^ 



F = . 



Introducing a new variable: 



one may rewrite (4.21) as follows: 

'v 



\xp 
ih 



iE 



-Vd 2 v -d v --)F(r,) = 



Finally, defining L(rj): 
one finds: 



F( V ) = e-^ 2 L(n) , 



L(n) = o 



(4.21) 
(4.22) 

(4.23) 

(4.24) 
(4.25) 



(iE 1 
V^7 " 2 

which is defining equation for Laguerre polynomials (cf. (3.9)). Hence we may define 



F+ = C n e-^ 2 L n ( V ) , n = 0,1,2,... , 



(4.26) 



where nth polynomial L n is given by (3.10) and C n is a normalization constant. The 
above formula for F+ is an analog of (3.11) for the oscillator Wigner functions. More- 
over, it follows from (4.25) that the spectrum is given by: 



(4.27) 



E n = ihj [ n + - ) , n = 0, 1,2, ... . 



Now, it is easy to check that functions F n defined by: 



F n :— Frt , 



satisfies 



H d * F n — F n -kH d — —E n F n 



(4.28) 
(4.29) 
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It follows immediately from the following property: 



f*9 = 9*f, (4.30) 

which may be easily proved using the definition of the ^-product (2.1). 

Now, let us study the basic properties of eigen-functions and compare these 
with those of oscillator Wigner functions W n . Clearly, contrary to W n are not real. 
Observe, that taking a constant C n in (4.26) according to 

(_iy» 

i.e. like in (3.11), one may prove that 

J F+(x,p) dxdp = J F~(x,p) dxdp = 1 , (4.32) 

in perfect analogy to W n . Moreover, F^ give rise to the following marginal probability 
distributions: 

J F±(x,p)dx = 5(p) , (4.33) 
F±(x,p)dp = 5{x) . (4.34) 



/ 



This property seems to violate the Heisenberg uncertainty principle - the particle is 
localized both in x and p variables. Clearly, we lose the probabilistic interpretation of 
F^ since the corresponding eigenvectors V'n (4-13) do not belong to the Hilbert space 
L 2 (R). Interestingly, F^ satisfy the following condition: 

Ft*Ft = ^ nm Ft , (4.35) 

in perfect analogy to (2.21). Therefore, one obtains the corresponding resolution of 
identity 

=2 E Rei? n = ^ ■ (4-36) 

n n 

Finally, it would be interesting to find relation between the resonant states 
defined in (4.13) and F^- Using some simple algebraic manipulations it is easy to show 
that 

F+(x,p) = C n j dye-™tf(x - \y)^~{x + \y) , (4.37) 

and 

F-(x,p) = Fj(x,p)=C n j dye-^^(x-^y)i^(x + ^y) 

= C n Jdye-™^(x-ly)il,+ (x + ly), (4.38) 

with C n defined in (4.31). Hence, each F^ is built out of ip^ an d VvT ■ Note, that 
these eigenvectors correspond to E n and E n , respectively. Clearly, if ip is a proper 
eigenvectors corresponding to a real eigenvalue E, then using the above prescription 
for F one recovers the Wigner function corresponding to if). Resonant states comes 
always in pairs and two members of each pair are needed to construct F. 
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Let us observe, that stationary functions may be denned in a more transparent 
way. Rewriting the Hamiltonian (4.9) as follows: 

H d (x,p) = -^(x*p + p*x) , (4.39) 

let us define Fq to be a normalized function satisfying the following conditions: 

p * F+ = and F$*x = 0, (4.40) 

which are solved by: 

F+(x,p) = \ e~ 2ixp / h . (4.41) 
Having the "+ ground state" F + one defines "+ excited states" by: 

F+(x,p) oc x n *F+(x,p)*p ri . (4.42) 
Analogously, let us define F ~ to be a "— ground state" satisfying: 

x * Fq = and F *p = 0. (4.43) 

One finds 

F (x,p) = — e 2 ^' h = F+(x,p) . (4.44) 

TTtl 

The corresponding "— excited states" read: 

F-(x,p) oc p n *Fv(x,p)*x n . (4.45) 
Using canonical commutation relation 

x*p — p-kx = ih, (4.46) 
one easily finds that F^ do satisfy: 

H d *F± = F±*H d = ±ih 7 (n+^J . (4.47) 

5. Harmonic oscillator vs. damped system 

Comparing the spectra of harmonic oscillator and damped system considered in the 
previous section one finds striking similarity, that is, they are related by the following 
relation: 

lo = ±Z7 . (5-1) 
Note, that performing the following canonical transformation: 

x = 75 {x + p) and p = 71 {x ~ Ph (5 ' 2) 

one obtains 

H d = - 1X p=^{P 2 -X 2 ) , (5.3) 

i.e. in the new variables (X, P), H d corresponds formally to the harmonic oscillator with 
co = ±27. This correspondence may be easily seen by observing that both Hamiltonians, 
i.e. (P 2 + X 2 ) and (P 2 — X 2 ) are related by the following ^-exponential: 

A A 2 

V x :=Exp(XXP/h) = 1 + -XP + ^XP*XP + ... , (5.4) 
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with A£l. Indeed, one may show that 

V x * X * V-x = e' iX X and V x * P * V- X = e iX P . (5.5) 
The above formulae imply: 

V A * (P 2 - X 2 ) * V^ x = e 2iX {P 2 - e~ 4iX X 2 ) , (5.6) 
and hence, for A = ±7r/4, one obtains: 



V ±tt/4 



1(P 2 -X 2 ) *V Tw/i = ±^(P 2 + X 2 ) , (5.7) 



i.e. both systems are related by a complex scaling V±^/4- Therefore, it should be clear 
that the corresponding eigen-functions W n and are also related by V± n /4. Let us 
denote 

2 

Now, if W n is an oscillator Wigner function satisfying: 



Ht = ±^{P 2 + X 2 ). (5.8) 



H± *W n = W n * H± = ±E n W n , (5.9) 
with E n given by (4.27), then F± defined by: 

F± :=V^ /4 *W n *V± n/4 , (5.10) 
satisfy the corresponding eigen-problem for the damped system: 

H d * F± = F± * H d = ±E n F± . (5.11) 
Moreover, it follows from (5.10) that 

F- = F£, (5.12) 

provided W n is real, and 

J dxdp = J (V± 7V /4 * V r T7r / 4 )W n dxdp = J W n dxdp , (5.13) 

i.e. F^ are normalized on R 2 . 

Finally, let us observe that introducing on R 2 polar coordinates (r, ip) the corre- 
sponding Koopman operator reads: 

L H ho = ivdp , (5.14) 

and hence, oscillator Wigner functions are SO (2) invariant since L>H ho W n = 0. On the 
other hand using hyperbolic coordinates (s,x)', P = scoshx, X = ssinhx, we obtain 

L H d = ild x , (5.15) 

and hence F^ are SO (1, 1) invariant. Moreover, this observation implies that the 
corresponding ^-exponential U*(t) has for the damped system following form: 

EM-m)tH d ) = cQsh 1 (t/2) exp [-2(i/n)tanh(t/2)fl- d ] , (5.16) 

which follows from (3.20). Note that 

{H d ,F} M = ih{H d ,F} , (5.17) 

and hence, like for the harmonic oscillator, quantum and classical evolution are given 
by the same formulae. 



11 



6. Damped harmonic oscillator 

Consider now a damped harmonic oscillator described by the following equation of 
motion: 

x + 272 + kx = . (6-1) 

As is well known this system plays a prominent role in various branches of physics, 
especially in quantum optics. The above 2nd order equation may be rewritten as a 
dynamical system on M 2 

±1 = -7x1 + ojx 2 , (6.2) 
X2 = —7x2 — wxi , (6.3) 

with ijj = \J k — 7 2 . Clearly this system is not Hamiltonian if 7 7^ 0. However, applying 
the procedure of [14] one arrives at the following Hamiltonian system on M 4 : 

x\ = {xi,H} = — 7x1 +L0X2 , (6.4) 

X2 = {x2, H] = -ujxi - 7x2 , (6.5) 

pi = {pi,H} = +7pi +ujp 2 , (6.6) 

p2 = {pi,H} = -ujpx + 7P2 , (6.7) 

where the corresponding Hamiltonian function is given by: 

Hdho{x,p) = u{pix 2 -P2X1) - j(piX! +P2X2) ■ (6.8) 

Let us observe that the above Hamiltonian may be rewritten as follows: 

Hdho(x,p) = <^(Pl *X 2 - P2*X\) - -(pi -kXi + X\ -kpi +P2 * X 2 + X 2 * P2) ■ (6.9) 

Now, let us introduce a new set of variables: 

X1+IX2 * Xx -IX 2 /c1nA 

ai = ' ai = ' (fU0) 



02 = -w a2 = ^r' (6 - n) 



satisfying the following commutation relations: 

{01, a 2 } M = {ai, a*} M = {02, aj^M = , (6.12) 

{ai,a2}M = {a2,a*}M = 1 ■ (6.13) 
Hamiltonian (6.9) takes in new variables the following form: 

H<±h.o = fi> {oi «2 * 0,1 + a * CL2 + w) 

= fta ( * ai + - J + fta ( * a2 + - ) , (6-14) 



where 

a = w — 27 . (6.15) 
Now, we are going to find the spectrum together with the corresponding eigenfunctions: 

H dho *F = F*H dho = EF . (6.16) 
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Define to be functions corresponding to "± ground states" , that is, 



°1 * F 00 = ) F 00 * a *2 = » 



and 



a 2 * F 00 = , F 00 * a\ = . 
Unique normalized solutions of (6.17)-(6.18) are given by: 

1 



and 

Moreover, denning 
and 

one shows 

and 

with 



^ 00 " (2nhy 6 



F oo - F oo ■ 



K m oc {<4) n *F+*a?, 

F nm « K) n *^00*< . 
Hdho * ^nm -^ram * -^dho F nm,F nm , 



Hdho * ^nm ^iitii * -^dho F nm,F nrn , 



(6.17) 
(6.18) 

(6.19) 

(6.20) 
(6.21) 
(6.22) 
(6.23) 
(6.24) 



i? nm = fta ( m + - ] — ha ( n + - ] = fkj(m — n) — ihq(n + m + 1) . (6.25) 



Let us compare the above formulation with the standard operator approach (see [1], 
[4]) based on the following Hamilton operator: 



#dho = oo{pix 2 - P2X1) - -{pix\ + x\p\ + P2X2 + X2P2) • 



(6.26) 



Following (6.10)-(6.11) we introduce (a^, at) which satisfy (6.12) and (6.13) with Moyal 
bracket { , }m replaced by the commutator. Now, let us introduce "± ground states" 
ip^ as the states satisfying: 



and 

— — ^* — n 
<W = a 2V 3 00 = U • 

Moreover, define two families of excited states: 
and 



It is easy to show that 
and 



<P nm :=o?(ai) m <^oo 



^dho fnm ~ F nm ( Pnm > 



Hdho fnm — F nmfr, 



(6.27) 
(6.28) 

(6.29) 

(6.30) 
(6.31) 
(6.32) 
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with E nm defined in (6.25). Using standard (xi, X2)-representation, i.e. x^ip = 
and pk = —ihdk^f, one easily solves (6.27) and (6.28). Up to non-important constants 
one obtains: 

<Poo(xi,X2) = S(x 1 )S(x 2 ) , (6.33) 

and 

^oo( x i' x 2) = 1 • (6.34) 

Clearly, neither nor ip^ belong to L 2 (R 2 ). Note, that there is a striking similarity 
between F^ m , ip^ m and F^, tp^ from section 4. Now, using a pair of resonant states 
(fnm one may easily show that 



Ft m oc J dy e - p >J m (x - |y)^ m (x + |y) . (6.35) 



Moreover, using a straightforward algebra one may prove the following 
Proposition. Normalized functions F^ m satisfy: 



<»Jx) = J F ± (x, p) dp = «5(x) , (6.36) 

^nm(P) = / F± (*> P) dx = 5 (P) > ( 6 - 37 ) 



and hence 



Moreover, 



I 



F± m (x,p)dxdp = l . (6.38) 



* f h - (2nh) 2 5nk ^ ml F ™ m ' (6-39) 
in analogy to (2.21). This implies the following resolution of identity: 

in accordance to (2.22). 

7. Another representation 

Both oscillator Wigner functions W n and the corresponding F^ and F^ m from sections 
4. and 6. respectively, are stationary function, i.e. they commute with the corresponding 
Hamiltonian. In the case of W n and F^ this property follows from that fact that 

W n = W n (H ho ) and F± = F ± (H d ) . (7.1) 

Now, in the case of a damped harmonic oscillator the corresponding Hamiltonian (6.8) 
may be written as a sum 

-f^dho = Hi + H2 , (7.2) 

where 

Hi = uj p A x and H2 = — 7P • x . (7.3) 

Clearly, 

F nm ^iim^) j (7 -4) 

and the stationarity of follows from 

{ffi,ff 2 } M = 0. (7.5) 
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Now, we show that it is possible to construct another family G nm such that 

G nm = G nm (Hi) , (7.6) 
and G nm satisfy the corresponding eigen-problem 

Hdho * G n m = G n m * -f^dho = l^nmG 'nm • (7 -7) 

Let us define Goo "ground state" satisfying 

a k *G 00 = and G 00 * 4 = , A; = 1,2. (7.8) 
Solving (7.8) one finds: 

G 00 = ef ^^-^^ = el xAp . (7.9) 

Clearly, Goo, contrary to F^, is not integrable over 1R 4 . Using (6.12) and (6.13) it is 
easy to show that the following set of functions: 

G nm = (air * (a* 2 ) n * Goo * a n 2 * a? , (7.10) 

satisfy 

Hdho * G n m = G nm * -Hdho = l^nrnG 'nm i (7-H) 

with 

Unm = ha + -^J + ha + ^ = hu>(n + m + 1) — ih^(n — m) . (7.12) 
Note, that 



G n rn — G mn , (7.13) 

and 

fJ'nm = H-mn • (7-14) 

Therefore, we have a natural pairing (G nm ,G mn ) in analogy to (F+ m , F~ m ). Interest- 
ingly, both approaches give completely different spectra of H&\ 10 : E nm and ^ n m defined 
in (6.25) and (7.12), respectively. 

Let us compare the above formulation with the standard operator approach (see [1], 
[4]) based on the Hamilton operator (6.26). The commutation relations may be easily 
represented in the space of functions of two variables (xi,j>2) : 

Sl = TTh { Xl - h w)> "* = Wn i Xl + h w) ' (7 - 15) 
22 = TTh (- p2 + h ^l) ' 25 = Wh i~ P2 " R £d ■ (7 ' 16) 



Introducing a ground state t/w 

aiV»oo = 02^00 = , (7.17) 

one finds: 

^o(xi,p 2 ) = e x ^ h . (7.18) 

Defining 

Y>nm = {a\) n {a* 2 ) m ^ , (7.19) 

one shows 

-f^dhoV'nm = /^rarnVVim i (7.20) 

with fi nm given by (7.12). One may show that it is possible to construct G nm defined in 
(7.10) out of resonant states Y> nm - However, contrary to F^ m , G nm are not normalizable 
and the striking analogy with Wigner functions is lost. 
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8. Concluding remarks 



In the present paper we analyzed the quantization of simple classical damped systems: a 
toy model defined by x = —jx and the damped harmonic oscillator. Both systems give 
rise to resonant states and the corresponding energy spectra are discrete and complex. 
It turns out that resonant states appear always in pairs: if ipi corresponds to E then 
there exists ^2 corresponding to E. We showed that each pair of such states may be 
used to construct an analog of the stationary Wigner function. Actually one constructs 
a pair of stationary functions 

dye~ %vv ^ x (x — -y) ^{x + -y) and F . 

A slightly different approach to quantization of damped oscillator was applied in [15]. In 
the forthcoming paper we show that both approaches are closely related. In a different 
context a quantum damped harmonic oscillator was recently analyzed in [16]. 
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